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Abstract 



We consider a number of uniqueness questions for several wide classes of active scalar equa- 
tions, unifying and generalizing the techniques of several authors [3JJ GH3 IH1 HI 131 HI] ■ As special 
fN) | cases of our results, we provide a significantly simplified proof to the known uniqueness result 

for the 2D Euler equations in L 1 n BM O [32] and provide a mild improvement to the recent 
results of Rusin [27] for the 2D inviscid surface quasi-geostrophic (SQG) equations, which are 
p I ■ now to our knowledge, the best results known for this model. We also obtain what are (to our 

\ knowledge) the strongest known uniqueness results for the Patlak-Keller-Segel models [SI HI IB]- 

We obtain these results via technical refinements of energy methods which are well-known in the 
'■^ | L 2 setting but arc less well-known in the H^ 1 setting. The H^ 1 method can be considered a 

generalization of Yudovich's classical method (34J [35] [25] and is naturally applied to equations 
such as the Patlak-Keller-Segel models with nonlinear diffusion, and other variants [4] [6] |3j . An 
important tool in our analysis is a Sobolev embedding lemma which shows that velocity hclds v 
with Vv G BM O are locally log-Lipschitz. 

> 

cn ! 1 Introduction 

The term active scalar refers to a wide variety of nonlinear advection-diffusion equations for a scalar 
density which, through a potentially nonlocal operator such as an elliptic PDE, gives rise to the 
advective velocity field. The purpose of this paper is study the uniqueness of sufficiently regular 
weak solutions to such equations in M d for dimensions d > 2. Our main interest will be achieving 
uniqueness at relatively low regularity and we will not be overly concerned with how quickly our 
solutions decay at infinity. 

The best known active scalar equations are the 2D Euler equations [y = 0) and Navier-Stokes 



equations (y > 0) in vorticity form 



uj-t + v ■ Vo; = !/Aw, 
v = V ± (A)- 1 w. 



The classical proof of uniqueness due to Yudovich [34] proves that weak solutions are unique provided 
w(t) G L oo (0,T;L 1 (M 2 )) n L°°(0,T; L°°(R 2 )). He also developed the classical existence theory for 
solutions to 2D Euler of this form based on the method characteristics (see also [25J ) . Given two weak 
solutions u)i,u)2 his method estimated the evolution of ||t>i(i) — f 2 (*) H2 (with the obvious notation). 
Measuring \\coi — W2||jj-i, however, provides a natural generalization of Yudovich's method to many 
other active scalars of interest (see [261 [6] 0] [3] and below). Later, Yudovich extended his proof to 
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cover cases with unbounded vorticity [35], but with certain restrictions on the growth of as 
p — > oo which ruled out logarithmic singularities. In that article, Yudovich also seems to have been 
the first to note the connection between the Osgood condition for the uniqueness of ODEs and the 
uniqueness of the 2D Euler equations. Vishik [32] later proved what is to our knowledge the optimal 
known uniqueness results which include the class u(t) € L°°(0, T; L 1 (IR 2 )) n L°°(0, T; BMO(R 2 )) 
(and other logarithmic improvements on L°°). His methods, in which of course the Osgood condition 
plays a role, are based on wavelet decompositions and a precise Littlewood-Paley decomposition 
of solutions to the 2D Euler equations in standard Eulerian form. One of the several results of 
our article is that Yudovich's significantly simpler method can also be used to prove uniqueness 
to the 2D Euler equations with u € L 1 n BMO. Ultimately, two relatively simple facts make 
this possible: that Calderon-Zygmund-type singular integral operators map BMO to itself |30| and 
vector fields v with Vv 6 BMO are locally log-Lipschitz (Lemma [T] below), and hence satisfy the 
Osgood condition. 

The first class of active scalars we study are PDE of the form, 

p t + V-( P Vp) = AA(p), (1) 

where A is non-decreasing and V is a linear operator which is roughly smoothing of order one (e.g. 
derivatives of solutions to second order elliptic PDEs). As mentioned above, the most famous of 
such equations are of course the 2D Euler and Navier-Stokes equations. Other examples include 
the parabolic-elliptic Patlak-Keller-Segel model and related variants [9l [3] where one would 
have for example V p = —V(A)~ 1 p, as well as the inviscid aggregation equations [22j 0]. As 
alluded to above, we will essentially apply Yudovich's energy method to these equations, measuring 
the difference of weak solutions in H . This also has the very useful additional benefit of treating 
nonlinear filtration-type diffusion very naturally, which is well-known to be a monotone operator 
in H -1 [31]. Indeed, the proof applies regardless of how pathological or strongly degenerate the 
nonlinear diffusion is. 

A second class of active scalar we consider are essentially variants of the 2D surface quasi- 
geostrophic equations studied by, for example [HI [331 QUI [2Ql [191 I27J , 

p t + Vp ■ Vp = -K-A) 7 P, (2) 

with v > 0, 7 € (0, 1] and V is a linear operator of roughly Calderon-Zygmund type which is 
divergence free, i.e. V • Vp = 0. The 2D SQG equations are the special case given by Vp = 
V J -(— A) -1 / 2 p = (— R2P1 Rip), where Rj are the standard Riesz transforms. To study these equa- 
tions we measure the difference of weak solutions in L 2 , but largely we use similar methods and 
techniques as used for the class ([TJ. Uniqueness results for both ([I]) and (|2j) are discussed below in 

m 

For inviscid active scalar equations, the authors are not aware of any uniqueness results which 
extend to problems with velocity fields less regular than those that satisfy the Osgood condition. 
Quantitative stability estimates, such as those here and those in [24} |23"1 [271 132] . transfer the 
regularity of the velocity field into stability in time of the norm being measured. In the Eulerian 
framework of an H~ l estimate here, or in the work of Vishik [32], the intuition for this is not 
immediately apparent. However, the Lagrangian methods of Loeper |24 [ I23 [ [TT ], [5]. which measure 
stability of inviscid active scalar problems in the Euclidean Wasserstein distance, make direct use of 
the Osgood condition as one which ensures that characteristics are well-posed and hence the induced 
flow map has sufficient stability and regularity properties. In contrast, non-uniqueness results for 
the 2D Euler equations [281 [29l HH] only apply to problems of far lower regularity than those which 
satisfy the Osgood condition, leaving a rather wide gap between the cases which are settled. The 
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DiPerna-Lions theory of renormalized solutions provides uniqueness for linear transport equations 
under very weak hypotheses, however, perhaps not coincidentally, these solutions do not seem to 
have well understood regularity and stability properties [T5| H]. 

For viscous problems, it seems fairly well-known that linear diffusion should improve the picture 
by providing a potentially useful negative term in stability estimates. In $4] we show that this can 
be applied also in the H~ l method for models of the form ([]]) with A(p) = p. As a special case, our 
results show that J * f \p\ 2 + \Vp\ 2 dxdt < oo is sufficient for uniqueness of weak solutions to ([I]). 

Some of the simplicity of our method arises from the use of scale-invariant (homogeneous) norms: 
for the first class we estimate stability in H , for the second class, we estimate stability in 1? . The 
use of more norms which could treat specific length-scales with more precision, such as Besov-type 
norms as in |32] , likely allows one to loosen requirements on the integrability of solutions at infinity 
and also treat active scalars such that the velocity operator V is of a type in between those above, 
for example as studied in [13]. All of this, of course, would complicate the proof significantly. 

Turning back to the 2D Euler equations as the canonical example, we point out here that 
the regularity implied by w £ BMO, as opposed to the integrability, is very crucial. Lemma [2] 
below shows that ui G L 1 n BMO imparts a certain amount of control on the integrability, that 
is, ||w|| p < p as p — > oo. However, we emphasize that this integrability alone is not enough for 
the following proof to work. Yudovich's proof in [35] uses only integrability and requires a stronger 
restriction on how fast \\u)\\ p can be allowed to grow in p. The results in [18] indicate that initial data 
which satisfies Yudovich's restrictions (but are not BMO) can result in velocity fields which are not 
log-Lipschitz (but still Osgood), and hence are in fact less regular than those studied here. Hence, it 
is important to emphasize that simply having logarithmic singularities in the density is not the key 
idea. The regularity conveyed by BMO is used twice in our work: that Calderon-Zygmund-type 
singular integral operators map BMO to itself and that velocity fields with Vv 6 BMO are locally 
log-Lipschitz (Lemma [T]). 

Notation and Conventions 

We briefly mention the definitions of weak solution we are using. Our theorems will generally require 
stronger regularity assumptions than these definitions; these simply provide the baseline. We use 
the following definition of weak solution to ([1]). This notion is stronger than a standard distribution 
solution, as we must be able to measure the stability of solutions -ff -1 , which will require taking 
test functions in H 1 . In 2D, a slight modification (suggested in [4]) must be made as if 1 (IR 2 ) is not 
a sensical function space. 



Definition 1 (Weak Solution to CD)). If d > 3 let V = H 1 (hence V* = H' 1 ) and if d = 2 
then let V = {/ € L°° : V/ G L 2 }. We say a measurable function p : [0, T] x R d ->■ R is a weak 
solution on some time interval [0, T] to if p t E L 2 (0, T; V*(R d j), VA(p) G L 2 (0, T; L 2 (R d )), 
pVp G L 2 (0,T;L 2 (R d )) and for all G L 2 (0,T;V) we have, 



For ([2]) we may take the standard notion of distribution solution as we need only to measure 
stability in L 2 . 

Definition 2 (Weak Solution to ([T])). We say a measurable function p : [0, T] x R d — y R is a weak 
solution on some time interval [0,T] to © if p G L 2 (0,T;L 2 ) and for all G C~((0,T) x R d ) we 
have, 
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We denote the set {u > k} := {x € D : u{x) > k}, if 5 C K rf then |5| denotes the Lebesgue 
measure and I5 denotes the standard characteristic function. We denote Ja -=j- A fdx. We denote 

the standard Lebesgue spaces ||/|| p = (J \f\ p dx) l ^ P and the weak Lebesgue spaces L p '°° associated 
with the standard quasi-norm 

\ Vp 

supa p A/(a) I , 

where Aj(a) = |{/ > a}| is the distribution function of /. We take the unitary convention of the 
Fourier transform, 



p,oo 



and denote the homogeneous Sobolev spaces H s as the closure of the Schwartz space under the 
norm H/ H^p := || |CT / an d the inhomogeneous Sobolev spaces H s as the closure of the Schwartz 

space under the norm ||/||_f/s := ||(1 + |C| 2 )' s ^ 2 /IIl2. 
We use M to denote the Newtonian potential: 

log \x\ d = 2 

f(d/2+l) 1 ,2— d . ^ o 

d(d-2)ird/2 \ X ~ y\ a - 6 ' 

We will denote space-time norms with the short hand 

\v t Ll = ll/Hl?(0,T;£S 



as the time interval and domain will usually be obvious from context. By weighted Young's inequal- 
ity we mean for a, b > and 1 = - + ~ and e > 0, 

ab < e p h e q — . 

P Q 

In formulas we use the notation C(p, k, M, ..) to denote a generic constant, which may be different 
from line to line or even term to term in the same computation. In general, these constants will 
depend on more parameters than those listed, for instance those associated with the problem such 
as /C and the dimension but these dependencies are suppressed. We use the notation / < p .fc,... g to 
denote / < C(p,k, ..)g where again, dependencies that are not relevant are suppressed. 



2 Lemmas 



We begin with a few relatively short lemmas which, while simple, provide the key technical tools for 
this work. The lemmas are simple enough such that we cannot guarantee they do not already exist 
elsewhere in the literature, but the precise nature of their connection to the relevant PDE problems 
seems to have gone overlooked. 

The first lemma is essentially a logarithmic variant of Morrey's inequality, allowing one to assert 
the log-Lipschitz regularity of a vector field based on the gradient having bounded mean oscillation. 
This lemma reflects probably the most fundamental observation, as this regularity is ultimately 
what provides the stability necessary to prove the uniqueness theorems below. 

Lemma 1. Let v G L p (M. d ) for some 1 < p < 00 such that Vv S BMO(R. d ) and 

sup / |Vu| dz < 00. 
x m d Jb(x,i) 
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Then v € L°° and for \x — y\ sufficiently small (depending on sup,j, effi d f B / x ^ | Vu| dz but not on 

x,y), 

\v(x) - v{y)\ < \x - y\ |log \x - y\\ \\Vv\\ B mo- 



Proof. We assume v £ «S(R ) and use density to extend to more general functions. Let x, y € 
and let r = \x — y\. Define W := B(x,r) n B(y,r). By the triangle inequality, 

\v(x) — v(y)\ < f \v(x) — v(z)\ dz + + \v(y) — v(z)\ dz. 
J w J w 



Moreover, 

C f 

\v(x) - v(z)\dz < — / \v(x) - v(z)\ dz, 

w \B(x,ry 

where C depends on the ratio between \B{x, r)\ and \W\, which is fixed in r. A standard fundamental 
theorem of calculus argument (see e.g. [16] pg. 267) implies 

\v(x)-v(z)\dz< [ J^^dz.. (3) 

JB(x,r) \X — Z\ 



\B(x, r) 



We are concerned with controlling this integral for r C 1. Without loss of generality we may assume 
x = 0. For notational simplicity, define S/v(y) = f(y) and B^ := B(0,2~ k ). Now, 



i/wi & </ m±iM dz+fBM<[)r . 

JB(O.r) \z\ 



[ 

JB(0,r) \z\ JB(0,r) 



Let us now focus on the first term, as the latter term is uniformly 0(r) by the assumption on the 
local averages of Vv. By / £ BMO, for all k £ N we have, 



/ \f-fB k \dx<c2- dk \\f\\ BM0 . 

JB k 



Moreover, as / £ BMO, we have 

\fB k ~ fB k+1 \ <d \ \f\\BMO, 

which implies \ fB k — $Bi \ < c k\\f\\BMO- We come to the main estimate, which breaks the integral 
into successive length-scales, 

I/- /fill , f , . f |/-/b x | , 

°° k2~ kd 

<d ii/Ubmo ^2 2-fc(^-i) = \\f\\ BM ° ^2 k2 ~ k - 

k>\og 2 r k>log 2 r 

An elementary computation shows 

\ kx k 1 = — > x k = — - . 

~ — ' ax ~— ' ax 1 — x 

k>log 2 r k>\og 2 r 



This finally implies, 

—dz < || Vv||bmo(^ I log r I - r) 



B{x,r) \Z - X\ 
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for r sufficiently small depending on f B r x ^ |Vi>| dx, which is uniformly bounded. Hence, for 
sufficiently small, 

\v(x)\<f \v{x) - v(y)\dy + f \v(y)\dy 

J B(x,r) J B(x,r) 



< 



!(x,r) J B(x,r) 

7v ^dz + Cr^- d \\v\ 

in— 1 11 1 



LP(B(x,r)) 



IB(x,r) \Z - X\ 

< ||Vu|| B Mo(r |logr| - r) + Cr d ^l v ~ d \\v\\ p . 
This completes the lemma. □ 



The John-Nirenberg inequality [30] is one classical way to quantify the idea that functions in 
BMO have at most logarithmic singularities and in particular, asserts that any / G BMO is in 
Lf . The following interpolation- type lemma provides a slight variant of this idea which will prove 
to be a simple but very important technical tool. 

Lemma 2. Let f G L Po n BMO, po < oo. Then, for p, po < p < oo we have, 

ll/llp SdP\\f\\BMO + II/IIpo 
In particular, for large p, the LP norm grows at most linearly as p — >■ oo. 

Proof. We will use the Calderon-Zygmund decomposition which splits / into a sum of two parts: a 
"good" part that is uniformly bounded piece, which we will bound using the pQ-noim of /, and a 
"bad" part which, though not bounded, has mean zero locally and may be bounded by the BMO 
norm of / (see [30j 1.4]). 

Without loss of generality, we may assume / has compact support. Define the dyadic maximal 
function 

Mf(x) = snp^- [ \f\(y)dy, 

where the supremum is taken over all dyadic cubes Q containing x. 

Let a > to be chosen later, E a = {Mf > a}, f\ = fxE c , f2 = f — fi- Note by the Lebesgue 
differentiation theorem f\ < a. Then 

||/i||£<a p - po ||/||fS- (4) 
Let {Qk} be the collection of maximal dyadic cubes with disjoint interiors such that 



\f\>a. 



Each x G E a is certainly contained in such a cube Q by the definition of E a , and that cube must 
necessarily be contained in E a , for otherwise, there is y G Q Cl E^, and by definition of Mf(y), 
Mf{y) > a, a contradiction. This maximality also implies 



\f\< d a. 

Since / is in BMO, we also have (see [Ml IV.1.3 (13)]) 

\f-f Qk \ P <d(p\\f\\ 



BMO) 
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We now estimate 



[\m p = Y,[ i/r<2 p £(7 \f-f Qk \ p +[ \f Qk \ p ) 

J U JQk U \JQk JQk / 



< d 2^Y J (p\\f\\BMoY\Qk\ + ^\Q k \) 

k 

= 2P((p\\f\\ BMO y + an\{Mf>a}\) 



IPO 

< d 2 P Mf\\BMo) p + aV)^ , 
where in the last line we used the Chebichev and Hardy-Littlewood maximal inequalities in showing 



i { M/> a} |<H£K<Ji/ir 



Pa 



a P0 ~- qjPO 



We now let a = \ \f\\ Po and obtain 



W < d 2V{{p\\f\\ BM o) p + \\f\\ p ) (5) 

Note that (jl|) implies ||/i||p < ||/||j> - Adding this estimate to ([5]) and taking pth-roots of both sides 
gives the final result. □ 

The final lemma provides straightforward sufficient conditions to be in the Hardy space % , 
which arises when treating the M 2 in Theorem [T] below. The lemma is also general enough to find 
use elsewhere [2]. 

Lemma 3. Let f € L 1 n L p for some p > 1 and satisfy f fdx = 0, M\ = f \x\ \f(x) \ dx < oo. 
Then f € H 1 and 



Proof. We will prove this using the duality between BMO and T-L . Therefore, let ||/C||.bmo = 1. In 
what follows, we denote KLr := |^^ ^ f^ <R K{x)dx. By the mean-zero condition on /, 

f JCfdx < [ \fC-Ki\\f\dx+ [ \K-Ki\\f\dx. 

J J\x\<l J\x\>l 

:=T 1 +T 2 . 

By Holder and p > 1 with ||/C||sAfO = 1, we have 

Ti < \\f\\ P (j^ ^\1C - JCxf dx^J 

Sp' II/IIp- 

The second inequality can be found in the proof of the John-Nirenberg inequality in [30j . We now 
deal with the second term. Define the dyadic annuli A n := {x G U. d : 2 n < |x| < 2 n+1 } for n > 0. 
Let E n := {x € A n : \K — JC\ \ > 2 n }. By definition we have, then 

T 2 = E/ |£-£i||/|<k+ / \lC-K 1 \\f\dx = T 21 +T 22 
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The second term can be estimated via, 

T 22 <2 n f \f\dx< f \x\\f(x)\dx. (6) 

Since K, 6 BMO, we have |/C 2 «+i — fC 2 n\ 1 ? and therefore \JC 2 n —K\\ <d n - Applying this to T 2 \ 
implies, 

T21 < (^j E \K-ldf dx^j /P \\f\\ p 

/ roo \ 1/p' 

P ' I |{|/c-/Ci| > x}nA n \\p'- 1 dx) \\f\\ p 

i/p' 



< (p' J \{\K-K 2 n\ > X-Cn}nA n \X p '~ l dX^j 



p- 



for some C = C(d) > 0. Using the John-Nirenberg inequality [30], there is some c = c(d) > such 
that 



f \K-K x \\f\dx< (\A n \ f 



\ 1 /p' 

e - c{X -Cn) X pi~l dX 



< [ 2 nd e Ccn / e^ Ac A p '" 1 dA 



p- 



Therefore, for large n we have, 

[ \)C-lC 1 \\f\dx<2^ +p '- 1 ^ p 'e Ccn/p 'e- c2n/p '\\f\\ p . (7) 

J E n 

Summing ([6]) and (J7j) over n then implies 

T 2 < p \\f\\ p + Mx. 

□ 

3 Uniqueness of weak solutions to (CD) and ([2]) 

We will require that the nonlocal linear operator V satisfy the following three conditions. 
CI \\Vp\\ g < \\p\\p with l + | = ^i + |andl<p<g<oo. 

C2 The regularity requirement that W : BMO(W i ) — > BMO(M. d ) is a bounded linear operator 
which additionally satisfies ||Vy/|| p < for all p sufficiently large. 

C3 The H- 1 stability estimate ||V/|| 2 < H/llff-i- 

Remark 1. The first condition can be weakened to require that V : L qi — > L q2 for any 1 < qi,q 2 < 
00. 
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Theorem 1 (H 1 Method). Suppose V satisfies properties Cl-3 and consider the active scalar, 

Pt + V-(pVp) = AA(p), (8) 

where A : R — > [0, oo) is measurable and non- decreasing. Suppose there exists two weak solutions 
p 1 (t),p 2 (t) defined on [0,T] which satisfy Pi (t) G L 2 (0, T; BMO(M. d )) n C t ([0, T]; Lg° w/iere */ 
d > 3 we assume 1 < po < 2d/{d + 2) and if d = 2 we tafce po = 1- In d = 2 we additionally assume 
Pi {x)x G Lf(0,T;L l x (R d )). Then z/pi(0) = p 2 (0) tfien pi(t) = p 2 {t) on [0,T]. 

Proof Define w := p\ — pi and let <f> = —N* w. First, notice by Lemma [2] that pi(t) G L 2 (0, T; L?) 
for all p > po- If d > 3, note also that Po/(po — 1) > 2d/(d — 2). By Young's inequality we have, 



||V0(t)||oo < ||V7Vl Bl (o)||_2d_||w||2d+ ||V^/"l Rd v fl (0) ||_pq_ 

which implies Vcfi G L 2 L^°. If (i > 3 we similarly have, 

| < ||^'lfli(0)||_M_l|«'||2d+ l|VA/'l R d\ i j l(0) 11 



po ' 



PQ- 1 



and therefore 0(i) G L 2 L£°. In M 2 , by Lemma O G rl 1 uniformly for a.e. time, hence, 

<f>(t) G Lf>L™ from the duality of V. 1 and BMO, 



\<f>(t,x)\ 



2tt 



log \x - y\ w(y)dy 



In particular, —A(p defines a bounded distribution and hence we have —A(p = w in the weak sense. 
In d > 3 note by Young's inequality for L 3771 ' 00 , 



l|V0(t)|| : 



< 



Ml 2d iivjvii d 



d+2 



hence V(p G L 2 L 2 . Proving this in d = 2 is more delicate as it does not follow from only L p estimates 
on w; we follow a procedure similar to that taken in [4]. By Parseval's theorem it suffices to prove 
that V0 g L 2 L 2 . First, notice 



\V(f>(t,x)\ 2 dx 



v<Kt,0 



dc 



l£l 2 



Mt,OI 2 ^- 



The latter term is controlled by Parseval's theorem since w G L 2 L 2 . The first term is not a 
priori controlled as |£|~ 2 is not integrable in R 2 , however, it suffices to prove that w(t) is uni- 
formly Lipschitz continuous. By the Riemann-Lebesgue lemma, w(t,x) G Ct([0,T]; L X (R 2 )) implies 
w(t,£) G Cf([0, T]; C X (M 2 )). Since ([8} is in divergence form, we have f w(t,x)dx = and therefore 
w(t, 0) = 0. The Lipschitz continuity follows from the control on the first moment: let £1,^2 G 
then by the mean value theorem for i G {1,2}, 



1 



< 



(2^) d / 2 
1 



(e"^ 1 - e~ ix <^ dx 



(2vr) d / 2 

< l£i - 6 



r/.r 



•Mi(ft) 16-61 ■ 
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Hence, V € l? x uniformly in time. Therefore, regardless of dimension we have that V0 € 
L^L 2 . n L 2 L^ and we can finally compute the evolution. By the regularity assumptions of the 
solution and the regularity properties of 4> deduced above, we may use as a test function in the 
definition of weak solution and compute the time evolution of ||V</>(t)||2, 



J \^^\ 2< ^ x = j ^ 4> ' ^ 4>tdx = j (ftwtdx 

= J <j)V ■ (VA( Pl ) - VA{ P2 ) + wV P i + p 2 (V Pl - Vp 2 )) dx 

= - j (A( Pl ) - A{p 2 )) (pi - P2 )dx - j • wV Pl dx - y V<£ • p 2 {V Pl - Vp 2 )dx 



:= TI + T2 + T3. 
Since A is non-decreasing with have, 

T\ = - J (A( Pl ) - A{p 2 )) ( Pl - P2 )dx < 0. 

We now deal with T2, which is most related to the regularity of the advective velocity field. Since 
w = —Acp and the boundary terms vanish due to the integrability conditions on w, 

^ J di(f>dij<j){Vpi)j + di4>dj(j)di(Vpi)jdx 

ij 

J d i (I l^l 2 ) (Vpih + d^dj^diivp^dx 

\Wpi\ \V<p\ 2 dx. 
Using VVpi € L 2 BMO x and Lemma [2] we then have, 

j \VV Pl \ \V4>\ 2 dx < \\VV P1 \\ P (J IV^I 2 ^- 1 ) dx 

^(pllvypillBA/o + llvypiiuo^iv^e/^ 1 )- 2 ! |V0| 2 ^ 
^(pllvypiiiBA/o + ^llv^ii^iiv^n 2 - 2 ^. 

Now we turn to T3, which regards the stability of the map p 1-4 V p. By Lemma [2] again we have, 



< 



1-1/p 

l-i/p 



|T3| < ||p 2 ||p||V0-(y /Ol -y /92 )|| p/(p _ 1) 

< (p\\P2\\bMO + ||P2||p ) ||V0[| 2 p/(p-l)||Vpi - Vp 2 \\2p/(p-l) 

< (p\\ P 2\\bmo + 1) ||v0||^||vpi - u P2 ||^||V0||^ 1/p ||y Pl - v P2 t 2 xlv . 

Assumptions CI and C2 together with Morrey's inequality imply Vpi € LjL^°. Assumption C3 
implies ||U(/9i — p 2 )\\ 2 < ||V0||2- Therefore, for all p sufficiently large we have 

d 

di' 



;||V^)|| 2 <p/(t)||V^)|| 2 - 2/P , 
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for some positive f(t) G L 2 ^ 1+2 ^ p \ Integrating, 



l|V#)|| 

Since pi(0) = p 2 (0) this implies, 



2/p 
2 



\\vm\\l /p 



< 



f(s)ds. 



t y 
f(s)ds) 







£? /2 (o,T)* 



1/3 



r3 /2, n ™. we have l|V0(i)||| < 2 p we may send p 

L t (U,J ) 



oo and 



Hence, if we restrict t < to < 2 31 

prove ||V</>(i)|| 2 = a.e. in time G [0, to], however clearly we can simply iterate and prove that 
||V0(t)|| 2 = on [0,T] and therefore w(t) = 0. □ 

The following result combines the methods of the above theorem with an L 2 estimate to provide 
a simplified and more general proof of the results in [27]. In particular, we do not need to assume 
that logarithmic singularities in the derivatives occur on a bounded set. 

Theorem 2 [L 2 Estimates). Let V be a bounded linear operator V : BMO(M. d ; E) ->■ BMO(R d ;R d ) 
which satisfies ||V/|| P < p||/||p for all 2 < p < oo and V • V f = 0. Consider the active scalar. 



Pt + Vp- Vp 



for some 7 € (0,1] and v > 0. Suppose there exists two weak solutions pi(t),p2(t) defined on 
[0,T] */i satisfy Pi (t) G C t ([0, T]; L 2 X (R d )) and Vpi G L 2 (0, T; BMO(R d ) DL Po (R d )) for some 
Po G [l,oo). TTien if pi(0) = p 2 (0) t/ien = p 2 (*) on [0,T]. 

Remark 2. It will be evident from the proof that we need only assume Vp G Ll + (0,T;BMOnLP°). 

Remark 3. If v > then we may potentially relax the conditions of this theorem; see $2] below. 
However, note that depending on the decay of Vp, there are still dissipative cases covered by this 
theorem that are not covered by Theorem [3] below. 

Proof. We estimate the L 2 norm of w := pi — p%. First, notice that by LemmaHJ pi, w G L 2 L^ and 
are in fact log-Lipschitz continuous for a.e. time. Using the divergence free condition, 



J \w\ 2 dx = — v J 



— — I 

2dt 



< 



I 



w(—Aywdx - 
I |V| 7 w\ 2 dx + 
wVp2 • Vwdx. 



wV P2VW + wVp2 ■ Vwdx 
wV P2 ■ Vwdx 



Using Lemma [2] for p > po and the boundedness of the mapping V : L 2p ^ p ^ — > L 2p ^ p 1 \ 



lit I 



i-i/p 



wVw\ p- 1 dx 

£ (p||Vp 2 ||sMO + ||Vp2||po) 



_2p_ 



(p-l)/2p 



2p 



(^||Vp 2 ||ijAfo + ||Vp 2 |U 
(pI|V/9 2 ||bmo + ||Vp 2 



2p 
p — 1 



2p 



(p-l)/p 



(p-l)/2p 



2/p|| w ||2-2/p 
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Since V>j G Lf (0, T; BMO x n Z|°), for p sufficiently large we have, 

i>K<p/( t )Hir" 



for some < /(i) G L^ 1+2/p ^ (0, T). As above in the proof of Theorem [H this suffices. □ 



4 Stability Results for Active Scalars with Dissipation 

In Sj3l the sign-definite dissipative effects were ignored and viscous and inviscid systems were treated 
in the same fashion. It is well-known that for the dissipative 2D quasi-geostrophic equations one 
can still have uniqueness at lower regularities than what is provided by Theorem [2] (see for example 
\6'i\ [TO] and the references therein). In light of existing results on 2D SQG, the following easy 
general theorem is proved by a standard application of Gagliardo-Nirenberg or Sobolev embedding 
and essentially mirrors contraction mapping arguments for local existence and small data critical 
theory, but we state and prove it for completeness. 

Theorem 3 (L 2 Estimates for Dissipative Active Scalars). Suppose there is some q > d/{2^) with 

2q 2q 

q > 1 such that V : Li- 1 — > Li- 1 is a bounded linear operator. Consider the dissipative active 
scalar, 

pt + v P • Vp = -K-A)V, 

for some 7 € (0, 1] and v > 0. Suppose there exists two weak solutions pi(t),p2(t) defined on [0, T] 
which satisfy Pi {t) G C t ([0, T]; L 2 (R d )) nL 2 (0, T; H 1 ) andV Pi G L r t (0,T; L q x ) with r = Then, 
if q > d/(2 r y) we have the stability estimate, 



\\pi(t) - P2(t)h < exp \ v *n-*C(d,q, ||Vpi|| L r £ «, ||^||i2 9 /(g-i)^ L 2,/( 9 -i))J ||pi(0) - p 2 (0) || 2 - (9) 

If d > 27 and the critical norm satisfies \\Vpi\\ d < i , eo(d,\\V\\ L 2d/(d-2 1 )_ > . L 2d/(d-2-y)), for some 

absolute €q depending only on those quantities indicated, then the solution map is a contraction in 
L 2 , that is, for all t > s > 0, 

\\pi(t) - P2 (t)\\ < ||piO)-p 2 ( S )|| 2 , Vt >s. 

Proof. Beginning as before we define w := p\ — pi and using the divergence free condition, 
1 d 



2dt 



J \w\ 2 dx = —v J w^A^wdx — J wV P2VW + wV p 2 ■ Vwdx 
V\ 1 w\ 2 dx+ I wVp2-Vwdx, 



where |V| is the Fourier multiplier |V| / := (2ir) d l 2 J |£| f{i)e lx '^d^. First consider the subcritical 

_d_ 

2-, ■ 



case q > 4-. By Gagliardo-Nirenberg we have, 



1 d d 

2 97 II IV7I7 „„ll 2 97 



Ml-*s_ ^ C(rf,9,7)lkll 2 297 II I v l 7w ll2 

9-1 



which implies, 

1 d f 

-— / \W\ 2 dx < —U\\ |V| 7 W|| 2 + ||W|| 2q \\VW\\ 2q IIV/02IL 

2 dt J q-l 9-1 



< — v\\ |V| 7 W\\2 + \\w\\ 2 2q WW 2? 29 II V/32 ||c 

< -HI |V| 7 w||l + C(d J g J 7)||«>||r"|| |V| 7 to||JlF|| 29 29 HVpall,. 

Li-' 1 -s-Lq-i 
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By weighted Young's inequality, 



7;^ [ H 2 dx < w \\l+ (7;) ^ d C (d,q,1f, ll^ll la_ In} 

2 dt I 2 \2J V l<7-i^l9-i 



279 



<U 2 " q - d cU q , 7 ,\\V\\ j. j. ||Vp 2 |||-- d ||«;||I. 

which implies the assertion ([9]) by Gronwall's inequality. 

We now prove the critical regularity assertion ([3]). Beginning as before we have, 

-— J \w\ 2 dx = — v J w{—/S.ywdx — J wVp2Vw + wS/p2 • Vwdx 

|V| 7 w\ 2 dx + / wVp2-Vwdx 



< -^WlVl 1 W\\l+ \\w\\ 2 _m_\\V\\ _M_ 2d ||Vp : 



^27" X,d-27^.^d-27 27 



By (homogeneous) Sobolev embedding we have, 

-- / \w\ 2 dx = -v\\\X7\^ w\\ 2 2 + C( 7 ,d)\\V\\ 2d 2a\\Vp 2 U\\\Vf r w\\l 
2 at J L d - 2 ~> -^L d - 2 f 27 



Hence, provided, 



C(j,d)\\V\\ 2d m ||Vp 2 |U < 1/, 

^d-27 _>.£d-27 27 



the solution map is a contraction in L . □ 

Using the Sobolev inequalities to capitalize on the stabilizing effects of the diffusion was straight- 
forward in the previous theorem. A similar (but slightly less obvious) approach can also be adapted 
for the H~ 1 method used above in Theorem [TJ The proof is limited to linear diffusion since, due to 
the lack of an appropriate embedding theorem, there does not seem to be an obvious way to extract 
a stabilizing effect from the nonlinear diffusions. One way this theorem is useful is that it shows 
solutions to ([1]) with "finite energy dissipation," J * f \p\ 2 + |V/?| 2 dxdt < 00 are unique (see Remark 
H] below). However, this theorem is obviously not as powerful as the deep methods of |17| . which 
ultimately prove that measure-valued solutions to the 2D Navier-Stokes equations in vorticity form 
are unique. A primary aspect of these methods is the approximation of the contribution of the 
atomic part of the initial data as self-similar Oseen vortices, an approach which could prove difficult 
or impossible in general. 

Theorem 4 (H~ l Method for Dissipative Active Scalars). Suppose, if d > 3 that there is a q > d/2 

2q 2q 

or if d = 2 that there is a q > 1 such that W : L q —¥ L q and VV- : L'J- 1 — > Li- 1 are both bounded 
linear operators. Consider the viscous active scalar, 

Pt + V ■ (pVp) = vAp, v > 0. 

Suppose there exists two weak solutions pi{t),p2{t) defined on [0, T] which satisfy 

2d 

Pi (t) e L£(0, T; L q (R d )) n L 2 (0, T; L 2 ) n C t ([0, T]-Lt 2 (^ d )) 

with r = 2qj(2q — d). In d = 2 we further assume u(x)x € L£°(0, T; (IR 6 ')). Then we have the 
stability estimate 

\\ P l(t) - P2(t)\\H- 1 <exp [v-^C(d,q,V)\\pi\\lr Lim ) |M0)-M0)||^-i. (10) 
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If d > 2 and the critical norm satisfies \\pi(t)\\d/2 — ve o{d,V) for some £q > depending only on 
the quantities specified for a.e. t G (0, T), then the solution map is a contraction in H , that is, 
for all t > s > 0, 

\\Pl(t) ~ P2(t)\\s-l < \\pi(s) ~ p 2 (s)\\ ti -i. 

Remark 4. Many solutions one can construct will not only have pi G L 2 (0,T; L 2 ) but also have 
V pi G L 2 (0, T;L 2 ). In this case, Sobolev embedding (in d > 3) or Gagliardo-Nirenberg (in d = 2) 
imply that pi G U t (0, T; L^) for a range of admissible pairs (r, (/) and hence all such active scalars are 
unique with only the assumption of pi G L 2 (0, T; H l )f~\Ct([0, T]; £ 2rf /( rf + 2 ) (]]j rf )) (with the appropriate 
modification in 2D). An interesting side note is that in 2D, Theorem |4] is only slightly stronger than 
Theorem (TJ as H 1 (R 2 ) 5MO(R 2 ) n I 2 (R 2 ) [21]. However, sometimes one is interested in 
solutions which do not have finite energy dissipation, for example the properly dissipative weak 
solutions studied in [7] . Unfortunately, it does not immediately follow from Theorem U that such 
solutions are unique; they are known to satisfy p G L^ 2 L^ 2 , but not p G L\L?J 2 . 

Proof. We prove the assertion (fTUj) first. Let w := p\ — p 2 and — Acp = w, equivalent to —N*w = <p 
since as in the proof of Theorem Q] we have <fi G L£L^°. Moreover, V0 G L 2 L 2 and computing the 
evolution of ||V^>(t)||2 similar to Theorem [fl 

\~dt I ' V ^' 2 dx ~ ~ V / H 2 dx + 2 J \Wpi\ \V^\ 2 dx + J P2 V<j)-Vwdx 

< -v\\w\\% + 2\\VVp 1 \\ q \\V<j>\\ 2 j 3 _ + ||p2|| ff ||^H|jia_l|V0||_aa_ 

9-I 8-1 9-1 

< -i/||u;||^ + 2||Vy|| L ^ L9 ||pi|| g ||V^|| 2 _^ + H/92 || 9 ||^V- || 2 Q 2 , ||V</>||V 
By the Gagliardo-Nirenberg inequality and the Calderon-Zygmund inequality, 

< d wml'^WD^wt < iiv^nr^iiA^ii* = imil^Mt 

9-1 

Hence, 

\jt j l v ^| 2 ^<-^II^Hi + CC^, V^H^illqllV^llJ - ^ II^Ill + C(^, V^V-)||P2|U[|V0I|2 - ^ ll^lll- 



By weighted Young's inequality 

1 d f ,„ ,,9 , . „ ,, i 

2dt 

Hence, by Gronwall's inequality. 



j \V0\ 2 dx < ~\\w\\ 2 + v~^C(d,q,VV,VV-) (\\ Pl \\i q9 - d + 11^11!^) ||V^|| 2 . 



|V^)|| 2 < ||V0(O)|| 2 exp 



u-^C(d,q,VV,VV-) f {\\pi(sW q +\\p 2 (s)\Qd S 

Jo 



This proves the first (subcritical) assertion. 

The second (critical) assertion follows similarly, with modifications to deal with the criticality. 
Beginning as before, 

1 d 



2dt 



J \V<j)\ 2 dx<-v j \w\ 2 dx + 2 J \S7Vpi\\V(f)\ 2 dx + j p 2 V<\>-Vwdx 
\ \w\ 2 dx + 2\\W\\ Ld /2^ L d/2\\pi\\ d / 2 II V^H 2 ^ 

J d-2 



< 

+ ||P2||d/2ll^V • || L 2d/(d-2)^ L 2d/(d-2) ||V^|| 2 2d 

d-2 
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By Sobolev embedding and the Calderon-Zygmund inequality again, 



ld_ 

Ydt 



\^<\>fdx < -i/||u>||f + C(d) ( ||W|| d a || d + || VV 




L(d-2) _>.£,T3=2j 



2d 2d 



2d 1 1 p2 1 1 d 



) 




Hence, the result follows immediately. 



□ 
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